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Nonlinear dynamical systems possessing an invariant subspace in the phase space and chaotic or 
stochastic motion within the subspace often display on-off intermittency close to the threshold of 
stability of the subspace. In a class of symmetric systems, the intermittency is symmetry-breaking 
[Ying-Cheng Lai, Phys. Rev. E 53. R4267 (1996)]. We report interesting and practically important 
universal behavior of robustness of supersensitivity, resonance and information gain in this class 
of systems when subjected to a weak modulation. While intermittent loss of synchronization may 
be harmful to application of high-quality synchronization of coupled chaotic systems, the features 
reported here may lead to interesting application of on-off intermittency. 



PACS number(s): 05.40.-a, 05.45.-a 



I. INTRODUCTION 

Nonlinear dynamical systems possessing an invariant 
subspace are of great interest, particularly when the 
system motion within the subspace can be chaotic or 
stochastic. Examples include chaotic systems with sym- 
metry or coupled chaotic systems || . Bubbling 
and on-off intermittency |^| are typical behaviors in sys- 
tem close to a threshold of transverse stability of the sub- 
space due to the fluctuative nature of the local transverse 
Lyapunov exponent in different part of the subspace. In 
bubbling, the subspace is stable on average, but local in- 
stability can result in large bursts away from the subspace 
when it is perturbed. The sensitivity of this weak stabil- 
ity to parameter mismatch and noise has been studied by 
Pikovsky and Grassberger (?]]. Intermittent loss of syn- 
chronization H in experiments with inevitable noise and 
parameter mismatch is undesirable in application of high- 
quality synchronization, such as to communication ||. 
In on-off intermittency, the subspace is slightly unstable 
on average, but local transverse attraction may keep the 
dynamics very close to the subspace for a long period of 
time. Great attention has been paid to the duration of 
this laminar period which exhibits universal power law 
distribution in a broad class of systems |(3) . Noise in the 
system prevents its state from approaching the subspace 
close beyond the noise level, thus has important effects on 
the laminar period distribution Q or the escape time [^0) . 
In noisy environment, bubbling and on-off intermittency 
are essentially the same phenomenon. 

Sensitivity in nonlinear systems can be very useful for 
applications such as controlling global dynamics of the 
system by local tiny perturbations Q. It is interesting 
to ask whether the sensitivity of on-off intermittency may 
lead to any potential application of the phenomenon. An 
observation is that in a class of symmetric systems, on-off 
intermittency can be symmetry-breaking jlj , namely, the 
bursting behavior does not possess the system symmetry 



when the system has two symmetric but distinct attrac- 
tors. This Letter reports that a combination of the sensi- 
tivity of on-off intermittency and the symmetry-breaking 
of the bursting can result in remarkable features in the 
systems subjected to a weak modulation signal in the 
noisy environment. 



II. THE MODEL 

Let x(t) represent the distance of the dynamics from 
the invariant subspace, and x(t) > and x(t) < denote 
the dynamics in the two symmetric components respec- 
tively. In the symmetry-breaking systems, transitions be- 
tween x > and x < can only occur when x(t) comes to 
the level of the weak noisy signal. For system displaying 
appreciable laminar state, main features of the dynamics 
can be described by the general linear equation close to 
the subspace: 



x{t) = [A + a^(t)]x(t) + cr 2 e(t) + ps(t). 



(1) 



Here A is the transverse Lyapunov exponent of the sub- 
space, and <Ji£(i) with (£(t)) = is the fluctuation of the 
local Lyapunov exponent due to the chaotic or stochastic 
motion within the subspace. In general, chaotic system 
has quickly decaying correlation, and in a large enough 
time scale t, has an asymptotic Gaussian distribu- 
tion. e(t) is the additive white noise with level o-i <C <J\ 
and s(t) is a weak modulation signal. The exact form 
of the signal is unimportant for the phenomena reported 
below, provided it varies on a time scale slower than the 
characteristic times of the systems. Here we consider s(t) 
a random binary stream (±1 with probability 0.5) with 
a bit duration T. p « ai is the amplitude of the sig- 
nal, and R = p/<J2 provides a natural measure of the 
signal-to- noise ratio (SNR). 

In Ref. Cenys and Lustfeld studied the statisti- 
cal properties of escape time of on-off intermittency sub- 
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jected to noise by means of Fokker-Planck equation. It 
has been shown that on-off intermittency is very sensi- 
tive to noise P,^0|, We employ the same approach of 
Fokker-Planck equation, and focus on the property of 
amplification of the weak external signal s(t) in the sys- 
tem. Our results will demonstrate that the system is also 
very sensitive to the weak signal, and the amplification 
of the weak signal is robust to the additive noise. This 
sensitivity exhibits resonant behavior as the system pa- 
rameters change. 

The Fokker-Planck equation for Eq. ([!]) is 



dW 
~~dt 



_d_ 

dx 
I d 2 
2dx 2 



(\ + -±)x+ps(t) 



[(a 2 x 2 + a 2 )W] 
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(2) 



In general, it is quite difficult to solve Eq. (g) exactly. 
Let To be the relaxation time of the system after s(t) 
switching from +1 (-1) to -1 (+1). If T > T , the 
probability distribution W(t) can establish an approx- 
imate static state during each bit of the input signal. 
Under the adiabatic approximation T ^ To, dW/dt w 0, 
and the static solution can be obtained analytically as 



W(x) = C [x 2 + -f 



2\ («-l)/2 



exp 



2ps(t) <j\x 
arctan 



(3) 



where a = 2\ja\. For \(Tix/A\ > 1 (A = max(p, a 2 )), 



irps(t) 



-sgnx 



(4) 



Now we see that the behavior of the system can be 
divided into two regimes. One is \x\ 3> |02e(t) + ps(t)\, 
where the dynamics is governed approximately by x(t) = 
[A + <Ji£(t)]x(t). Let z = In \x\, then z(t) = A + <n£(t) 
which describes a Brownian motion with a constant drift 
A and diffusion constant o\j2. The nonlinearity of the 
system can be modeled by an effective reflecting bound- 
aries of the Brownian motion at ix b , which is of the 
order of <j\. The probability density has a power form 



W{x) 



IQ-1 



but is asymmetric for x > and x < 



in the presence of s(t). The system can rarely perform 
transition between x > and x < in this regime due 
to the symmetry-breaking property, until it comes to the 
other regime, where the noisy input a^eit) + ps(t) dom- 
inates the dynamics and the system performs transition 
between x > and x < frequently. The behavior of the 
system is determined by the competition between the dif- 
fusion and the drift of the Brownian motion. If the drift 
time tj, = ln(<7i/A)/|A| is much smaller than the diffusion 
time td = 2m 2 (ai/A)/a 2 , the drift is dominant over the 
diffusion, and the system will either come to a metastable 
state induced by the noisy input for A < 0, or approach 
some state away from the invariant subspace for A > 0. 



In both cases, the weak noisy input has no significant ef- 
fects on the system behavior, i.e. the system is insensitive 
to the modulation. On the other hand, the diffusion is 
dominant for t b 3> id, and the system can have access to 
both the level of the weak input and the boundary of the 
nonlinearity, exhibiting typical on-off intermittency and 
sensitivity to the weak modulation. Variation of the pa- 
rameter A or tj\ affects the competition between the drift 
and the diffusion, and the system is expected to display 
optimal response to the weak modulation with resonant 
characterization. 



III. ROBUSTNESS OF SUPERSENSITIVITY TO 
THE WEAK SIGNAL 

Let us consider the diffusion dominant region close to 
the critical point of the stability of the subspace, e.g. 
|A| < 1, \aln(x b /A)\ < 1, 01 ~ 1, p ~ 10" m (m > 1). 
Employing the approximation in Eq. (||), and the effec- 
tive reflecting boundaries at ixb, the ensemble average 
(x(t)} is estimated as 



(*(*)} « «(*) 



x b itR 
tann , 



ln(x b /A) 



(5) 



which in the noise-free limit oi — > 0, assumes the form 



(*(*)) ~ «(*) 



ln(xb) — lnp 



(6) 



A logarithmical dependence of (x(t)} on the input level p 
means an amplification of the weak signal ps(t) with an 
factor (x)/p ~ 10 m /m (m ^ 1), i.e. the system exhibits 
supersensitivity to extremely weak modulation close to 
the critical point. This sensitivity was also reported in an 
overdamped Kramers oscillator with multiplicative noise 
free from additive noise (02 = 0), which is a specific 
example in this class of systems In the absence 

of s(t), the system produces symmetric bursting pattern 
with (x(t)) = 0; while the bursting pattern is reorganized 
to manifest the weak signal after it is fed into the system 
(see Fig. 1). The most interesting and practically impor- 
tant property is that the weak signal is manifested even 
buried in a relatively high level of noise, namely, the ro- 
bustness of the supersensitivity. This behavior originates 
from the symmetry-breaking of the on-off intermittency 
in the system. 

To demonstrate the above analysis, we employ the fol- 
lowing system in simulations |l|] 

y = -iv + %(i - y 2 ) + fo sinwt, 

x = (a + by) sin(x) — x + oie(t) + ps(t), (7) 

where y constitutes the forced Duffing chaotic oscillator. 
With 7 = 0.05, f = 2.3 and w = 3.5, the Duffing sys- 
tem is chaotic and o\ 1=3 0.9646. The nonlinearity of the 
variable x is related to an experimental model of super- 
conducting quantum interference device (SQUID) [O. 
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However, we should stress that the specific form of the 
nonlinearity is of no importance for the phenomena. The 
transverse Lyapunov exponent of the invariant subspace 
x = Q is X = a — 1 due to the symmetry of the Duff- 
ing chaotic attractor. Fig. 1 shows typical behavior of 
the system and good agreement between the analytical 
and the simulation results for (x(t)} as a function of R. 
The agreement demonstrates that the general stochas- 
tic model in Eq. (Q) gives good account for this type of 
system even though the motion in the subspace is deter- 
ministic chaos. 



IV. RESONANT BEHAVIOR 

It is difficult to perform generally a quantitative anal- 
ysis of the system response to the weak signal based on 
the linear dynamics in Eq. (Q) as parameters A or o\ 
changes, because the effective boundary Xb changes with 
the nonlinearity and the linear dynamical model with 
an effective reflecting boundary is often not sufficient to 
capture the dynamical property if a is appreciately posi- 
tive. Moreover, as the parameters change, the relaxation 
time T may become comparable to the bit duration T, 
and the transient behavior plays an important role in the 
system response and an adiabatic approximation is not 
valid any more. To demonstrate the resonant properties, 
we rely on simulations with the system in Eq. (Q), while 
the Brownian motion model can provide a qualitative 
understanding of the properties, thus showing that the 
properties are generic and universal for a general class of 
the systems. 

For a system with on-off intermittent output x(t), the 
ensemble average (x(t)) and the correlation between s(t) 
and x(t) is relatively small even for the noise-free case 
(72 = 0, due to the power law fluctuation of x(t). To bet- 
ter characterize the response of the system to the mod- 
ulation, we transfer the output series x(t) into a binary 
stream X(t) by a threshold crossing process: suppose x(t) 
becomes larger than a prescribed threshold Xth at some 
moment, after that X(t) will keep at X(t) = 1 until x(t) 
crosses — Xth at another moment; X(t) will not switch 
back from X(t) = — I to +1 until x(t) crosses Xth again, 
and so on. This binary presentation captures the most 
important feature of the transition of the bursting pat- 
tern between the two symmetric attractors. X(t) has a 
strong correlation with s(t) for weak noise case ai < P if 
the system is close to the critical point. The exact value 
of x t h is not crucial for the properties described below. 
In the following, we fix p = 10~ 7 , T = 2000 and Xth = 1, 
and take the cross-correlation function C between s(t) 
and X(t) estimated using 10 4 bits of a random stream of 
s(t) to demonstrate the resonant behavior in the system. 

(1) With the change of A. For A rather below the crit- 
ical point A = 0, the system has a metastable state close 
to the level of the noisy input, and the diffusion is not 
strong enough to produce large bursts frequently enough, 



C will be small. On the opposite, if A is rather above 
the critical point, the drift is also dominant so that the 
system can seldom access to the level of the weak mod- 
ulation, and becomes insensitive to the switching of s(t) 
between ±1, resulting in a small C again. Close to the 
critical point, the system can access to the level of the 
weak signal and produce large bursts frequently due to 
strong enough diffusion. The switching of the weak sig- 
nal is "sensed" and manifested by asymmetrical bursting 
to x > and x < 0, giving an optimal value of C. This 
precess is illustrated by C as a function of A in Fig. 2 for 
various R values. 

(2) With the change of <j\. For a small o\ where the 
drift is dominant over diffusion, the system is not sensi- 
tive to the weak noisy input and C assumes a small value. 
With the increase of 01, the diffusion becomes stronger, 
resulting in a smaller relaxation time To and more fre- 
quent large bursts, and the system becomes more sensi- 
tive to the weak input. In the noise- free case 02 = 0, 
the increased sensitivity enables X(t) to keep closer in 
phase to the weak signal s(t) and C approaches closer 
to 1.0 if the system maintains to work in the symmetry- 
breaking regime, and in general a resonant behavior is not 
expected. The picture becomes quite different if 01 7^ 0. 
With smaller To and increased sensitivity, the system can 
keep up with and manifest more and more noise-induced 
transitions in shorter time-scales, and the transition rate 
of X(t) between ±1 may become much higher than that 
of s(t), leading to a decreasing C. An optimal response is 
achieved when the diffusion is strong enough to become 
sensitive to the weak input but not too strong to mani- 
fest a lot of noise-induced transitions in short time-scales. 
Typical example of the system response as a function of 
(7 1 is shown in Fig. 3. 

The above resonant behavior is similar to the con- 
ventional stochastic resonance where a dynamical sys- 
tem displays increased sensitivity to a subthreshold sig- 
nal with an optimal level of additive noise, see Ref. Gfl 
for an extensive review. Resonance occurs when a noise- 
controlled time-scale in the system matches that of the 
signal. In our system, the underlying mechanism of the 
resonant behavior is quite different. The sensitivity to 
an extremely weak signal is induced by the multiplicative 
chaotic or stochastic motion in the subspace. To achieve 
this sensitivity, it is required that the system is in the 
on-off intermittency regime so that it can become suscep- 
tible to the weak signal by coming close enough to the 
subspace, and manifest and amplify it by quick enough 
large bursts away from the subspace with symmetry- 
breaking. As system parameters A and <7\ change, a 
competition between these two factors leads to the res- 
onant behavior. More interestingly, resonant behavior 
with respect to the change of ax, the level of the mul- 
tiplicative chaos (noise), necessarily occurs only in the 
presence of the additive noise due to the nature of this 
competition. These features are generic and universal in 
a general class of systems displaying on-off intermittency 
with symmetry breaking. Multiplicative stochastic res- 
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onance has been studied by Gammaitoni et al |15j in a 
multiplicatively driven bistable system with A = 1. In 
that case, the system is out of the regime of on-off inter- 
mittency, and consequently cannot display the property 
of (super)sensitivity, and the resonance with respect to 
the change of A was not resported. 



V. INFORMATION GAIN 



Now consider the system from the viewpoint of trans- 
mission and amplification of a weak signal ps(t) contam- 
inated with channel noise <72e(i) through a system dis- 
playing on-off intermittency with symmetry-breaking. It 
is very interesting and practically important that more 
information about the signal may be obtained from the 
output X(t) than from the noisy input ps(t) + <r 2 e(t) it- 
self, besides the fact that the weak signal has been ampli- 
fyed to a level discernible with a low resolution detector. 
To examine the information gain, we compare C with the 
correlation between the signal ps(t) and the total noisy 
input ps{t) + (T2e(t), i.e. Ci„ = R/Vl + R 2 . C, Ci n 
and their difference are shown in Fig. 4. C comes to a 
saturated value for R > 1, where the noise- induced tran- 
sitions between x > and x < in short time-scales 
are rarely manifested by large bursts. This value does 
not approach 1.0 as Cj n due to an average time-delay 
between X(t) and s(t) induced by the relaxation time 
which is longer for weaker signal. Clearly, appreciable 
information gain is obtained by the system in a broad 
range of i?, and an optimal gain is found at a certain R 
value. If taking into account the effect of the time-delay 
between X(t) and s(i), e.g. by defining C T as the max- 
imum of the correlation between X(t) and s(t — r), the 
information gain region can be wider (Fig. 4). 

To conclude, we demonstrate interesting universal fea- 
tures of robustness of supersensitivity, resonance and in- 
formation gain in a class of nonlinear system subjected to 
a weak modulation. These systems present a new mech- 
anism of resonant behavior compared to conventional 
stochastic resonance. While intermittent loss of synchro- 
nization may be harmful for any applications employ- 
ing high-quality synchronization Q, the features found 
in this letter are meaningful for potential applications 
of on-off intermittency. On-off intermittency has been 
demonstrated in many experimental systems and we be- 
lieve that the behaviors reported in this work can be 
tested in physical experiments. 
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FIG. 1. (a) A typical time series of x(t) in the noisy en- 
vironment with 02 = 10 -6 : upper panel, R = 0, and lower 
panel, R = 0.5. The dotted line is s(t). A = -0.02, T = 2000, 
(b) Numerical estimated (a;) as a function of R (stars) com- 
pared to the analytical estimation of Eq. (^) (line). A = 0, 
b = 1 and p = 10" 7 . 



FIG. 2. Resonant behavior with respective to A = a — 1; 
6=1. 

FIG. 3. Resonant behavior with respective to b. A reso- 
nance occurs when 02 7^ 0. 



FIG. 4. An illustration of information gain. X = 0, 6 = 1. 
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